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Abstract 

In this work, a minisuperspace model for the projectable Hofava-Lifshitz (HL) gravity 
without the detailed balance condition is investigated. The Wheeler-deWitt equation is 
derived and its solutions are studied and discussed for some particular cases where, due 
to HL gravity, there is a "potential barrier" nearby a = 0. For a vanishing cosmological 
constant, it is found a normalizable wave function of the universe. When the cosmo- 
logical constant is non-vanishing, the WKB method is used to obtain solutions for the 
wave function of the universe. Using the Hamilton-Jacobi equation, one discusses how 
the transition from quantum to classical regime occurs and, for the case of a positive 
cosmological constant, the scale factor is shown to grow exponentially, hence recovering 
the GR behaviour for the late universe. 
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1 Introduction 



Hofava-Lifshitz (HL) gravity is a quite original proposal for a ultraviolet (UV) completion of 
General Relativity (GR) p], in which gravity turns out to be power-countable renormalizable 
at the UV fixed point. GR is supposed to be recovered at the infra-red (IR) fixed point, as the 
theory goes from high-energy scales to low-energy scales. In order to obtain a renormalizable 
gravity theory one abandons Lorentz symmetry at high-energies [H [2] . Even though the idea 
that the Lorentz symmetry is a low-energy symmetry has been previously considered [3], the 
novelty of the HL proposal is that the breaking of Lorentz symmetry occurs the very way 
as in some condensed matter models (cf. Ref. pQ and references therein), that is through 
an anisotropic scaling between space and time, namely r — > bf and t — > b z t, b being a scale 
parameter. The dynamical critical exponent z is chosen in order to ensure that the gravitational 
coupling constant is dimensionless, which makes possible a renormalizable interaction. As the 
Lorentz symmetry is recovered at the IR fixed point, z flows to z — 1 in this limit. 

The anisotropy between space and time leads rather naturally to the well known 3 + 1 
Arnowitt-Deser-Misner (ADM) splitting pE], originally devised to express GR in a Hamiltonian 
formulation. Following Ref. p.], a foliation, parametrized by a global time t, is introduced. Since 
the global diffeomorphism is not valid anymore, one imposes a weaker form of this symmetry, 
the so-called foliation-preserving diffeomorphism. Choosing this approach, the lapse ADM 
function, N, is constrained to be function only of the time coordinate, i.e. N = N(t). This 
assumption satisfies the projectability condition pQ. In order to match GR, one could also choose 
N = N(f,t), a model dubbed non-projectable and which has been investigated in Refs. [5J |6]. 
The next step involves getting a gravitational Lagrangian into this anisotropic scenario. For 
this purpose, the effective field theory (EFT) formalism is used: every term that is marginal or 
relevant at the UV fixed point (z ^ 1) is included and, at the IR fixed point, only z = 1 terms 
survives. GR is then presumably recovered. The number of terms that must be included splits 
HL gravity into two classes, depending on whether one adopts the detailed balance condition 
or not. It is argued in Ref. pQ, that if one allows every relevant term to be included into the 
Lagrangian, the number of coupling constants would be so large that any analysis would become 
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impracticable. The detailed-balance condition is inspired by non-equilibrium thermodynamics 
[7] and, loosely speaking, it states that the potential terms of a D- dimensional action are 
obtained using a (D — l)-dimensional function, the superpotential. It is argued that although 
detailed balance is a simplifying assumption, it is by no means a necessary one [HI OJ - It is shown 
that the list of allowed terms is not so large after all and the detailed balance Lagrangian is 
obtained after a proper choice of coefficients. 

A common problem plaguing all HL versions is the presence of a scalar field mode, which 
has a trans-Bogoliubov dispersion relation with k 6 term JTUJ fTTj . This scalar mode appears, 
likewise a Goldstone boson, after the breaking a continuous symmetry. To avoid this mode, 
one has to introduce more symmetries: besides the foliation-preserving diffeomorphism, a local 
U(l) symmetry can be introduced [12] and it is shown that the scalar mode is then gauged 
away. This version of HL gravity is referred to as general covariant, given that the number of 
degrees of freedom matches the one of GR. For reviews on these versions of HL gravity, the 
reader is referred to Refs. [H [101 E] • 

Cosmological considerations have been extensively studied in the context of HL gravity (for 
reviews, see Refs. [131 Ej)- One subtle point that arises in the projectable version concerns 
the lapse function which being just a function of time, implies that the classical Hamiltonian 
constraint of GR is no longer local, and must be integrated over spatial coordinates. It is shown 
in Ref. [15] that this yields an additional term that mimics dust into the Friedman equations. 
However, as noted in Ref. [9], the Robert son- Walker metric is homogeneous, so this spatial 
integral is simply the spatial volume of the space and hence this "dark dust" constant must 
vanish [16] . The presence of terms higher spatial curvature terms in HL gravity, gives rise to a 
plethora of new cosmologies, that exhibit, for instance, some bouncing and oscillating solutions 
[T7l [TBI [19] . A classification of Friedman-Lemaitre- Robertson- Walker (FLRW) cosmologies 
for HL gravity has been performed in Refs. [T6l [20] . One should notice that the analysis 
carried out in these references is entirely classical and based on the resulting Friedman-like and 
Rauchaudhury-like equations. 

Quantum cosmology is an interesting step towards the understanding of quantum gravity 
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and the initial conditions of the universe. Its setup consists in splitting space-time, using the 
ADM formalism and applying well known quantum mechanical considerations for constrained 
systems. The cosmological principle is evoked so that space-time is foliated in leaves with 
a constant global time. To implement the quantum scenario, one promotes the Hamiltonian 
constraint H = to an operator equation, the Wheeler-deWitt (WdW) equation, Hip = 
0, where ip is the wave function of the universe [21]. The WdW equation is an hyperbolic 
equation on the space of all 3-metrics, the so-called superspace. Its complexity makes the task 
of obtaining solutions a formidable one. To deal with this equation, one often considers simpler 
spaces, such as for instance the FRLW space-times, which leads to a minisuperspace model, 
where the number of degrees of freedom is considerably reduced from infinite (any 3-metrics) 
to one (the scale factor) [22]. Despite their relative simplicity the minisuperspace models are 
not completely free from problems. Indeed, one can point out, for instance, the fact that 
the wave function of the universe is not in many cases normalizable, which implies that the 
usual interpretation of quantum mechanics cannot be used. However, in the context of some 
particular models normalizable wave functions have been found and discussed [23l [241 [25]. For, 
comprehensive reviews, see e.g. Refs. [26l [27] . 

We argue that quantum cosmology allows for a valuable insight of HL gravity in the quantum 
context. In both approaches, one foliates the space-time in constant global time leaves, a 
procedure that automatically satisfy the projectability condition. But when adopting the QC 
formalism in HL gravity, one faces the problem of turning the Hamiltonian constraint into the 
WdW equation, since the Hamiltonian constraint in the HL gravity is not local. Nevertheless, 
choosing a FLRW metric minisuperspace model or, more generally, a spatially homogeneous 
cosmological metric, one can argue that the spatial integration yields a local constraint. Notice 
that another suitable feature of HL gravity, is that it does not introduce higher than first 
order time derivatives of the scale factor on the action, making the quantization procedure 
straightforward as a mixture between time and spatial derivatives and powers of momentum is 
not found. Indeed, the kinetic part of the Hamiltonian has the same structure as the one in the 
usual QC, and HL gravity introduces only higher spatial derivatives terms, which dominate on 
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the very small scales. Notice that the problem of high order time derivatives imposes severe 
obstacles for applying QC on high order derivatives gravities and string theory, but not for HL 
gravity One concludes then that the minisuperspace model can be naturally implemented in 
the HL gravity proposal. In the minisuperspace model, the WdW equation for the HL gravity 
was obtained in Ref. [28], however in there the interest was on the cosmological constant 
problem and the HL WdW equation was neither discussed nor solutions were presented. 

In this work, one investigates the projectable HL gravity without detailed balance in the 
context of the minisuperspace model of quantum cosmology for a FLRW universe without 
matter. This particular choice, despite being much simpler than the non-projectable version 
j5j E] and the general covariant approach [12] , exhibits the main features of the HL gravity and 
contains the detailed balance as a limiting case. A matter sector is not introduced, given that 
the main interest in the very early universe, where the HL new terms dominate and for the 
late universe, an epoch dominated by the cosmological constant. Moreover, the inclusion of the 
matter sector and how it is coupled to HL gravity remains an open question [TO] . 

This paper is organized as follows: in section |2j the minisuperspace HL is presented and the 
WdW-equation is obtained. In section [3j the solutions of the WdW are found and discussed. 
In section HJ the wave function of a HL is interpreted and an analysis of the Hamilton- Jacobi 
equation is performed. Finally, concluding remarks are presented in section [5j 



2 The Wheeler-DeWitt equation 



2.1 Metric 



One considers the RW (Robertson- Walker) metric with R x S 3 topology 



ds 2 = a 2 {-N(t) 2 dt 2 + a 2 'y ij dx i dx?) 




where i,j = 1,2,3, a 2 is a normalization constant, N(t) is the lapse function and 7^ is the 
metric of the unit 3-sphere. Its metric is given by 7^ = diag (jzrz, r 2 , r 2 sin 2 6) . 



The extrinsic curvature takes the form: 
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K l3 = -i- ( + V.iV, + V.A, ) , (2) 



2a A V ^ 

where A* is the ADM shift vector and V, denoes the 3-dimensional covariant derivative. As 
Ni = for RW-like spaces in study, 

K« = -~9a- (3) 

Taking the trace one gets 

K = K ij 9ij = — — -. (4) 
yi ° aNa K J 

The Ricci components of the 3-metric can also be obtained, as the foliation is a surface of 
maximum symmetry 



fl=-5-5- (6) 



2.2 Hofava-Lifshitz action 

The action for the projectable HL gravity without detailed balance is given by [51 IS]: 



ShL = it / d * xdtN ^9 i K » Kij ~ XR2 ~ 9oMl - giR - 92M- 2 R 2 - gzM-fRijRV- 

-g 4 M-*R 3 - g s M-*R {R'Jl 1 ,) - ge M~ x 4 R^ R\ R\ - g 7 M;*RV 2 R - ^ 8 M P " 4 V \R jk V> R jk } , 

(7) 

where gi are coupling constants, M P1 is the Planck mass and Vj denote covariant derivatives. 
The time coordinate can be rescaled in order to set g\ = —1, recovering GR value. One also 
defines the cosmo logical constant A as 2 A = goM 2 v An important feature of the IR limit is the 
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presence of the constant A on the kinetic part of the HL action. GR is recovered provided A — > 1 
(corresponding to the full diffeomorphism invariance), however A must be a running constant, 
so there is no reason or symmetry that a priori yields A = 1 GR value. Phenomenological 
bounds suggest however that the value of A is quite close to the GR value [TO] . 
Performing these redefinitions the HL action reads 



S HL = —^- J d 3 xdtN^ {KijRV - \K 2 + R - 2A - g 2 M m 2 R 2 - g z M~ 2 R tj R ij - 

g 4 M-*R 3 - g 5 M- A R {R l 3 R\) - g^M^R^R* - g 7 M~ 4 RV 2 R - g 8 M^V 'iR jk V ; R jk } . 

(8) 

2.3 Hofava-Lifshitz minisuper space action 

In order to consistently reduce the number of degrees of freedom when restricting the 3-metrics 
of the superspace to be isotropic and homogeneous, one can consider that the restriction is 
performed directly into the equations of motion, or through the substitution of the RW metric 
into the Lagrangian density and then obtain the equations of motion for the remaining degrees 
of freedom. In general, the physical content of these two ways are different, showing that the 
restriction cannot be done over the Lagrangian unless one properly solves the arising constraints. 
For the RW metric, without matter fields, these procedures are shown to lead to the same results 
[29J. Since the HL proposal introduces only an anisotropy between space and time, it does not 
alter the homogeneity of RW metric and hence the metric ([1]) can be substituted into Eq. (JS} 
yield the HL minisuperspace action 



M p 2 , x 2?r 2 x 3(3A - l)a 2 f / -a 2 a 6a 2Aa 2 a 3 

ohl — z / dtN 



2 J { N 2 3(3A-1) 3(3A-1) 

12 24 1 ' 

- M " X 3(3A - 1)^ X (3 * + - M " X 3(3A - X ( " 4 + 395 + *>} ' 

where the spatial integral J cPx^^y = 2tt 2 has been performed. A further simplification is 
obtained after choosing units so to satisfy a 2 x 6n 2 x (3 A — l)M|j = 1. The minisuperspace 



action then reads 



S HL = - [ dtN f ~ d2a < 2a AM p " 2 a 3 24n 2 (3g 2 + g 3 



2 J { N 2 (3A — 1) 18tt 2 (3A - l) 2 a 
2887T 4 (3A-1)(9^ 4 + 3g 5 +g e " U0) 



a 3 



Following Ref. [16] the dimensionless coupling constants are redefined as 



9c 



9a 



3A- 1' 
AM: 



18tt 2 (3A- l) 2 ' 
g r =247i 2 (3g 2 + g 3 ), 

g s =288vr 4 (3A-l)(9(74 + 3( ? 5 + ^)- 



(11) 



Notice that gc > 0, which stands for the curvature coupling constant. The sign of #a follows the 
sign of the cosmo logical constant. These two terms are already present in the minisuperspace 
GR model, but now they depend on A. The coupling constants g r and g s can be either positive 
or negative as their signal does not alter the stability of the HL gravity (cf. Ref. [16]). As 
discussed in Ref. [9] , physically, g r corresponds to the coupling constant for the term behaving 
as a radiaton and g s stands for the term behaving as "stiff" matter (p = p equation of state). 
The minisuperspace action is finally written as [9] 



S 



HL 



dt 



N 



v , +9ca - g\a -g r ^ 



(12) 



2.4 Horava-Lifshitz minisuperspace Hamiltonian and Wheeler-DeWitt 
equation 

The canonical conjugate momentum associated to a is obtained using Eq. (TT2"]) 
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The Hofava-Lifshitz minisuperspace Hamiltonian density is performed using Eqs. (1121) and 



H = U a a- C = \- (-III - g c a 2 + g A a 4 + g r + . (14) 



The next step to implement the quantum cosmology programme involves promoting the 
classical minisuperspace Hamiltonian into an operator on which the so-called wave function of 
the universe is applied to [2TJ [22] . 

This is subtle point in HL gravity since there is no global diffeomorphism, just a foliation- 
preserving diffeomorphism [Ij. This can be also seen as the lapse function no longer depends 
on the space-time variables, as in GR, but now it depends only on the global time N = N(t), 
as discussed in Sec. [U This implies that the Hamiltonian constraint is not local, however this 
problem can be circumvented for an homogeneous metric, like Eq. ([I]), as the integration over 
space can be performed as seen above. The canonical quantization is obtained by promoting 
the canonical conjugate momentum into an operator, i.e. U a i— > —i4r- Due to ambiguities in 
the operator ordering, one chooses = — -^-^ { aP ~tk) [22J . The resulting WdW equation is 
then obtained 



The choice of p does not modify the semiclassical analysis [30J, hence one chooses p = 0, 
and the WdW equation is written as 



- g c a 2 + gAa 4 + g r + 2$ \ V(a) = 0. (16) 



da 2 ' ' a 2 

This equation is similar to the one- dimensional Schrodinger equation for h = 1 and a particle 
with m = 1/2 with E = and potential 
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V(a) = g c a 2 - - SV 




a 



(17) 



2.5 Hofava-Lifshitz minisuper space potentials 

The WdW equation derived in the last section resembles an unidimensional Schrodinger equa- 
tion with potential given by Eq. (|T7|) . Classically the allowed regions are such that V(a) < 
since E = 0. A complete analysis of the phase structure of HL FLRW cosmologies was per- 
formed in Ref. |16j . 

Notice that the first two terms of Eq. ( fTTl) are the usual GR terms of the quantum cosmology 
analysis [22]. HL gravity introduces the last two terms which dominate for a <C 1, i.e. are 
relevant at short distances, presumably at the very early universe, where the GR description 
must be replaced by the quantum gravity one. At the very early universe, this potential is 
dominated by the term —g s /a 2 , implying that for g s < 0, this potential exhibits a "barrier" 
that might prevent space-time to get singular. The case g s > is not examined as it leads to a 
cosmology that cannot be suitably investigated using QC techniques. Notice that the detailed 
balance condition yields g s = 0. 

The choice g s < splits the discussion into three distinct scenarios, for positive, negative 
and vanishing cosmological constant. In what follows one studies the cases ensued by these 
choices for the coefficients Eqs. (TTTj) . 

2.5.1 A = case 

In this case the curvature term dominates at large distances. The universe oscillates between 
ai and a 2 . The potential, depicted in Figured! is written as 



V A=0 (a) = g c a 2 - g r 




a 



(18) 



2.5.2 A y£ case 



For large a, the potential is dominated by the cosmological term, and given by 
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V A=0 (a) 




a 



Figure 1: Potential for A = 0. 



Va^o(o) = gco 2 - Qko 4 - 9r 




a 



(19) 



The sign of g\ follows the sign of the cosmological constant A. For positive A the potential 
is depicted in Figure |2a] and the negative A case in Figure |2b] For a positive cosmological 
constant, one considers a potential that has three positive roots (ai, a<i and 03) hence there 
are two classically allowed regions for a\ < a < a2 and a 3 < a and a forbidden region where 
02 < a < 03. There is another possibility, discussed in Refs. [131 [16], m which the potential has 
only one real positive root, namely a%. The expression for the roots of Eq. ( Tl9|) . can be found 
in Ref. [31] . and will not be presented here as their expresions will not play any role in what 
follows. The potential can be factorized as 



By the same token, for a negative cosmological constant, one finds a similar behaviour 
already present in the A = case: classically, the universe oscillates between a% and 02. This 
implies that the potential Eq. ( fl9|) reads 



VA>o(a) 



a 



A I 2 2\ / 2 2\ / 2 2\ 

- [a — a 1 ) (a — a 2 J [a — a 3 J . 



(20) 



V A <o(a) 



a- 



- [a 2 — af) (a 2 — a 2 ) (a 2 + a 2 ) , 



(21) 



where a = ±iaj are the imaginary roots of this potential and a, is real. 
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V A>0 (a) 



V A<0 (a) 




a 



a 



(a) A > Potential 



(b) A < Potential 



Figure 2: Potentials for non- vanishing cosmo logical constant. 



3 Solutions of the Wheeler-DeWitt equation 

Having described the three types of potentials one encounters, the task is to solve the WdW 
equation (fl6l) . If A ^ 0, the cosmological constant term is quartic so the Eq. (fl6l) cannot be 
solved in a closed form, and the WKB approximation will be employed. 

3.1 Boundary Conditions 

To find suitable boundary conditions for the WdW equation one has to rely on physical assump- 
tions. The most discussed choices are the de Witt boundary condition [21], the "no-boundary" 
proposal [22l [26] and the tunnelling boundary condition [32J . 

The deWitt boundary condition [21] is the one in which the wave function of the universe is 
required to vanish wherever there is a classical singularity. It is inspired on quantum mechanics, 
and it is suitable for the cases under study here, as there is a potential barrier (bounce) for 
a < 1 yielding that the singularity is inside a classically forbidden region (cf. Ref. [23] ). The 
de Witt boundary condition is expressed, for FLRW models as: 



The "no-boundary" condition, of Hartle and Hawking [22 j . arises from obtained using the 
Euclidean path integral formalism. In that formalism, the ground state for the wave function 
of the universe is written as (cf. [33]) 



ipdw(a = 0) = 0. 



(22) 
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V>(a) = / [da]exp(-I), (23) 
Jc 

where C denotes that the integral is taken over compact manifolds, and / is the Euclidean 
version of the action defined in Sec. [2] such that the corresponding Euclidean action / = —iShl 
can be obtained from Eq. (j!2p using dr = iNdt and the N — 1 gauge: 



-iS 



HL 



dr 



da\ . 3 , 9r . 9s 

dr J a a A 



(24) 



where r is the Euclidean time. It is possible to evaluate ip(a) nearby a = [33j [34]. In this 
case, it can be proved that for r < 1 (close to a = 0), one has ^ = 1 [26J. Substituting these 
conditions into the Euclidean version of action Eq. (124j) . and integrating from to Ar, one 
finds 



dr 



-(1 + g c )r + £ A r 3 + 



5r , 5^ 



r 



(25) 



which is / — > — oo, for ^ r 7^ and g s 7^ yielding a divergent wave function. This shows that 
the "no-boundary" condition is not suitable to the problem under study. 

It is important to notice that the boundary condition ip(0) = does not mean that there 
is a quantum avoidance of the classical singularity given that it is a sufficient but not a nec- 
essary condition [2"3"| |2"4"| |2"5] . In the above references, some examples are given where ip — > 0, 
but J da\ip(a)\ 2 diverges, and conversely cases where j da\ip(a)\ 2 — > 0, but ip diverges. The 
conditions under which the classical singularity is removed or avoided by quantum mechanics 
are understood only in specific cases (cf. j23j |25] and references therein). 



3.2 WdW equation solution for a«l 

This region corresponds to the very early universe, where the HL terms dominate. This HL 
epoch is expected since any theory of quantum gravity is supposed to alter the GR description 
of the structure of the space-time at small distances. For this case, the WdW Eq. (TIB]) reads 
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which is an Euler equation whose solution is ip(a) = a s and S = | ± |yT — 4<? s . As the solution 
must be real and ip(0) = 0, one finds that the wave function for a -C 1 goes as 



^(a) ~ a^^ 1 ^ 7 , (27) 

which yields that g s < 1/4, a "quantum" bound for g s . Notice that this coefficient is uncon- 
strained by classical consideration [16]. Although the potential for g s > corresponds to an 
infinite well, this quantum bound gives rise to a mild singularity, which admits a well-defined 
mathematical treatment (cf. Ref. [35J). 



3.3 WdW equation solutions for a ^> 1 

This limit corresponds to the very late universe, which is dominated by the curvature and cos- 
mological constant terms. These terms are already present in the usual GR quantum cosmology 
setup [221 EE], reflecting the fact that GR behaviour is recovered at large distances. For A = 0, 
Eq. f|T6|) for a ^> 1 is given by 



^-g c a 2 }Ha) = 0, (2J 
which has the following asymptotic solution 



^(a)~ e -^ a2 . (29) 

Thus, as expected, the wave function has an exponential behaviour, since a ^> 1 corresponds 
to a classically forbidden region for the potential Eq. ffTH]) . 

For the positive cosmological constant case (#a > 0), the WdW Eq. flT6|) reads 



d 2 

— +g A a i \4>(a) = 0, (30) 
14 



whose asymptotic solution is given by a combination of Bessel and Neumann functions with 
v = 1/6 (cf. Eq. 8.491.7 of Ref. [36]). Since the limit a > 1 is being considered and the 
Neumann functions only diverge at a = 0, these two functions are admissible: 



V(a) = C lv ^Ji (-|V) + CtVENi (-^a 3 ) , (31) 

where C\ and Gi are constants. A further analysis show that the asymptotic expansions for 
Bessel and Neumann functions of any order [y) and large arguments (\z\ — > oo) is given by (cf. 
Eqs. 9.2.1 and 9.2.2 of Ref. [SI]): 



N »( z )*y— sin { z ~ 



(32) 

2 . / VTX 7T 

y ~ i 



Substituting these asymptotic expressions into Eq. f l3~Tj) one finds the asymptotic behaviour 
of the wave function for large a: 



where Cj = Ci^J for z = 1,2. Notice that this wave function is oscillatory, denoting that 
this is a classically allowed region and damped as |^(o)| 2 ~ a~ 2 - Not surprisingly, the same 
behaviour is found in Ref. [22] when the cosmological constant dominates the evolution of the 
universe, and the GR regime is recovered. This issue will be discussed in Sec. HI Interestingly, 
the WKB method yields the same asymptotic expression given by Eq. (I3T1) [251 135] . 
For A < =^ g A < 0, Eq. ([16]) is written as 

{^-(-^A)a 4 }^H = 0, (34) 

whose asymptotic solution is a combination of the modified Bessel functions, l v {z) and K v {z) 
(cf. 8.406 of Ref. [36J) of order v = 1/6. However, I v {z) grows exponentially as z — > oo (cf. 
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Eq. (I3"6"j) ). hence only K„(z) represents an acceptable solution for large a. The wave function 
in that limit is given by 



V>(a) ~ ^aKi h^^a 3 ) . (35) 

As above, asymptotic expansions for modified Bessel functions of any order (u) and large 
arguments (\z\ — > oo) are obtained using Eqs. 9.7.1 and 9.7.2 of Ref. 



e 

Uz) 



^ (36) 



One then gets that the wave function for the very late universe 

—\/— OA ^ 

V>(a) ~ -e^^ a . (37) 
a 

Thus, one concludes that for a ^> 1 the wave function is, as expected, strongly suppressed in 
this limit, given that this region is classically forbidden. 

3.4 WdW equation solution for A = 

If the cosmological constant vanishes the WdW Eq. ([TBI) read|j] 



d 

i o :j< '• ;ir > 

aa z gt 



— -g c a 2 + g r + ^}V(a) = 0. (38) 



After a change of variables, x = g^a, Eq. (1381) reads 



2 + ^ + f H*) = °- (39) 



2 , _9V , 

si' 

This equation can be exactly solved in terms of the associate Laguerre functions. Indeed 
the following differential equation (cf. Eq. 22.6.18 of Ref. [3~T] ) 



1 In a quantum mechanical context, this equation was solved in the Problem 4 of §36 of Ref. [37] . 
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has as solution y(x) 



= e 2 x' 




Laguerre functions. Comparing Eqs. (1391) and (j4"U|) . one finds that 



a 



9r 



■An + 2 + 2a 



(41) 



V9c 



ip(a) 



where N is a normalization constant to be obtained below. As g s < 0, a > 0, the wave 
function ip(0) is regular. Comparing with Eq. ( 12T1) . one also finds that the ratio g r j \fg~c must 
be quantized. This is not surprising given that one is solving the Schrodinger equation with 
E = 0, which for a bounded potential has a discrete spectrum. Thus, E = is an eigenvalue only 
for specific values of the coefficients and these values must be quantized. Another interesting 
feature of this solution is that it is normalizable. Indeed, using Eq. ( l4"Tj) and Eq. 8.980 of Ref. 

/. I 1/4 
r(n+o+i) • ^ ne com pl e te solution, 

for the A = WdW equation satisfying the ^(O) =0 boundary condition is given by 



This wave function behaves as Eq. f l27|) for a ^ 1 and as Eq. (1291) for a ^> 1. The Laguerre 
associated function L" is an rt-th order polynomial which yields that the wave function has n 
nodes. It is not difficult to verify that for a fixed g s , any g r obtained through the quantization 
condition Eq. (14T]) implies that the potential Eq. ( !T8|) has two positive real roots and so the 
universe is oscillating whatever value of n > is chosen. For large g r values, which implies that 
n is large, the very structure of the wave function shows that the universe is almost classical 
(this can be also seen by inspection of the probability density plots for large n as shown in Figure 
IHdl) . This result can be understood in terms of Bohr's correspondence principle, according to 




(42) 
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II 1 I 

(c) n = 2 (d) n = 10 

Figure 3: Probability density for the wave function of the universe for diferent values of n. The 
dashed line plot represents the potential Eq. (fl8l) . 

which the classical behaviour is obtained from the quantum one in the limit of large quantum 
numbers. The probability density distribution for the wave function of the universe, |?/>(a)| 2 , for 
some n values are plotted in Figure [31 One clearly sees that the solution is highly suppressed 
in the classically forbidden region, and it is oscillating with n nodes in the classically allowed 
region. Finally, it is straightforward to show that the singularity is avoided in this model given 
that the probability to find the universe at a = vanishes due to the HL gravity terms. 



3.5 WdW solution for A > 



If 9a 7^ 0, Eq. ffl6l) cannot be analitically solved. The behaviour of the wave function for large 
a and nearby the singularity a = were already discussed. For the intermediate regions where 
the curvature term starts to become relevant, after the HL epoch (very early universe), one has 
to rely on the WKB approximation, which for the classically allowed region is given by 



if? (a) 



1 

WW 11 



exp 



±i / y/\V(a)\dc 



ai 



(43) 



where ± symbol denotes that one must consider a combination of these two exponentials, V(a) 
is the potential Eq. ([20]) and ai is the classical turning point, i.e. V(ai) = 0. The following 
integral must be solved 



'i, J a\ a 



One uses that V(a) < for the classically allowed region and, hence |V(a)| = —V(a). This 
integral is valid for a\ < a < a2 < a%. Changing the variables to t = a 2 and rationalizing the 
square root, one finds 



This integral can be written as a sum of elliptic integrals [36j [38]. Using the formulas 
3.131.3, 3.132.2, 3.137.3 of Ref. [36] and the reduction formula 230.01 of Ref. [38] one finds 
after a rather long although straightforward computation 



y/\V(a)\da 



9a 
3 



sj(a 2 - a\)(al- a 2 )(a 2 3 - a 2 ) - (a\ + a\ + a 2 ) £(7, q)+ 



2d\a\ + a\a 2 + a 2 ,^ — af 



F(l,q)- 



3a| a 3 



n 



a? — al 



(46) 



where 7 = arcsin \/ ; 9 



\_ \ . F((p,k) is an elliptic integral of first kind, E((p,k) is 
an elliptic integral of second kind and II (tp, n, k) is the elliptic integral of the third kind (cf. 
[36J). Substituting Eqs. ([20]) and (146]) into Eq. ([43]) one finds the WKB wave function of the 
universe. This approximation is valid only if dV J? <C |V(a)| 3 ^ 2 [37] . 



For a positive cosmological constant, the classically forbidden region, ai < a 2 < a < a 3 , has 
the following WKB wave function 



il>(a) 



Ci 



|V(a)|V4 
and one has to solve 



exp 



V|y(a)|da 



+ 



Co 



\V{a)\V* 



exp 



V|V(a)|da 



(47) 
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Following the same steps as before, this integral can be written as 



2 tV(t- ;)(t-oi)(o§-t) 

which can be solved using the formulas 3.131.5, 3.132.4, 3.137.5 of [36] and the reduction 
formula 230.01 of [38]. One gets 



P y/Vtfjda =^ { J (a* - aj)(a 2 - al)(a§ - J) - ( °° + $ t^j <® U (.,p 2 ,p) 

Ja 2 [ V a 3 ~ a l 

( aj - a\a\ - a\a\ + a 2 2 aj \ 3a 2 2 {a 2 2 - a\) ( p 2 a\ \ \ 

(50) 



// 2 2\( 2 2\ / 2 2 

where k = arcsin \ ,1 2^2 i\ and p = \ °% °| . Substituting Eqs. ( 1201 and (1301 into Eq. 

V l a 3~ a 2;l a ~ a lJ V a 3 a l 

(j47"|) . one finds the WKB wave function for a x < a 2 < a < 03. 



3.6 WdW solution for A < 



One is interested in the WKB wave function for the classically allowed region a\ < a < a 2 . The 
WKB wave function is given by Eq. (|43p . The steps are the very ones of the above procedure, 
however, following the discussion of Sec. |2l one must consider that the smaller root (a 2 ) is 
negative and real (cf. Eq. (l2~Tj) ). Using the formulas 3.131.5, 3.132.4, 3.137.5 of Ref. [36J and 
the reduction formula 230.01 of Ref. [3H], one gets 
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f yWi^da =^ I J(al - a 2 )^ - a?) (a 2 + a 2 ) - ^ + n {n,p\p) + 



of + a 2 a 2 + a|g| - 5a|a| \ ^ _ 3a 2 .(a 2 + a 2 y }J , 




(51) 



/ (a 2 —a 2 )(a 2 —a 2 ) la 2 — a 2 

where k = arcsin < / ; 2 _ ^ and p = w § ■ The WKB wave equation is obtained after 



(a 2 -a 2 )(a 2 -a 2 ) 

inserting Eqs. ( 12TT) and ( 1511) into Eq. ( 1431) . This completes the WdW solutions for the potentials 
given by Eqs. ( !T8l) and ( Fl9l) . Any analysis of the WKB wave function in this regime, due to its 
complex expressions is somewhat difficult. 



4 Intepretation of the wave function 

To analyse the behaviour of the wave function, one must compute K? = KijK % \ the trace of the 
square of the extrinsic curvature Eq. ((2]) [33], [33] - If a wave function is oscillatory (exponential), 
K? has positive (negative) eigenvalues. Using Eqs. ([3]) and ( !T3|) . one gets 

9 d 2 

= —TIT?- (52) 

cHa 4 da 1 

Defining the auxiliary quantity W := and the asymptotic expression Eq. ( 12"7|) . for 

a <C 1, one obtains that IV < 0. When the HL gravity terms dominate the universe, the wave 
function is exponential, corresponding to an Euclidean geometry. 

For the very late universe, the behaviour of the wave function is very different depending 
on the value of the cosmological constant. If A = 0, the wave function is given by Eq. (12T?|) and 
it is easy to show that W < 0. Showing that the behaviour is exponential. If A > 0, the wave 
function is given by a combination of oscillatory functions Eq. ( |33|) . giving W > 0, meaning 
that the geometry is Lorentzian or classical. Finally, for negative values of the cosmological 
constant, Eq. ( |3~T|) yields W < 0. In the case studied here, the computations are quite simple, 
and it is not surprising to find this result since the nature of the wave function given by Eqs. 
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(|2"7j) . ( 129|) . (13"3"]) and (|3"7j) can be obtained directly by inspection. 

The semiclassical approximation implies that the configurations will oscillate about the 
classical solution [27]. In order to verify whether GR can be recovered for the low-energy limit 
in this approximation, one obtains the Hamilton- Jacobi equation from the WdW equation 
through the WKB method. The WdW equation (cf. Eq. ([IB]) ) can be written as 

^-V(a)^(a) = 0; (53) 

substituting a wave function of the form ip = Re [CV S ] , where C is a slowly varying amplitude 
and S is the phase, one obtains the Hamilton- Jacobi equation 

^V + V(o) = 0. (54) 
da J 

Notice that S is real provided V(a) < (classically allowed region), denoting that the wave 
function is oscillatory. If V(a) > (classically forbidden region), S is imaginary and the wave 
function has an exponential behaviour. Using Eq. (154]) . one finds that S = J ^/—V(a)da is the 
phase for the classically allowed region. This integral was discussed in Sec. El Applying n a to 
the wave function of the universe, ip = Re [Ce lS ] , one gets H a ip = —i^- = —i + ^C^f ) z lS ■ 
The WKB assumption, |^| < yields that U a = Using Eqs. (fl3l) and (131]). one is lead 
to 

a(t) a 

da—===, (55) 

a(0) y/-V[a) 

this equation relates the global time and the scale factor. One intends to investigate the 
regions where the positive cosmological constant dominates, so V(a) ~ —g^a 4 , substituting 
this asymptotic potential into Eq. ( 155]) . one gets that the time evolution for the scale factor 
when a 3> 1 is given by 

a (t) ~ e^K 1 . (56) 
This corresponds to a de Sitter expansion phase, a behaviour expected for the GR regime, which 
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is recovered when A = 1. Of course, this does not prove that GR is recovered as an IR fixed 
point of the HL gravity, but shows that a HL FLRW cosmology yields for a>la semiclassical 
solution that corresponds to the GR one. A possible approach to tackle the problem would 
involve considering a scaling A = A(a), and expect that A — > 1 for a 3> 1. Considerations of this 
nature were developed for cosmological models with scale-dependent Newtonian gravitational 
coupling (cf. Refs. [391 HQ]). 

5 Conclusions 

In this work, the quantum cosmology for the Hofava-Lifshitz gravity without matter is inves- 
tigated for the closed universe. In the minisuperspace model, the WdW equation is derived, 
and it is shown that the HL gravity introduces terms that are dominant for short distances, 
modifying the behaviour of GR on these scales. One chooses the configurations for which the 
HL gravity new terms act as a "potential barrier" close to the singularity, a = 0. 

The solutions for the WdW equation are obtained considering the deWitt boundary con- 
dition, Eq. (|22|) . which states that the wave function vanishes at the singularity. For a < 1, 
corresponding to the very early universe when the HL gravity terms dominate, the wave func- 
tion is an exponential, that is typical of classically forbidden region. A quantum bound for the 
coefficient g s is found. 

For the very late universe, when o>l, the curvature and the cosmological constant terms 
dominate and, one finds that, for A = or A < 0, the wave function is exponentially suppresed, 
denoting as before, that this region is classically forbidden. For a positive cosmological constant 
case, one finds a damped oscillatory behaviour, already found in the usual QC for GR. 

For the vanishing cosmological constant, an exact solution can be obtained. In this case, 
one finds that the singularity is avoided due to quantum effects as the probability to reach 
the singularity a = vanishes and that g r is quantized. Fixing the value of g s , for large 
values of n (large g r values), one can obtain a classical universe according to the analog of the 
correspondence principle of old quantum mechanics. The complete exact solution for A ^ 
cannot be obtained, although wave functions in the WKB approximation can be obtained for 
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the intermediate regions. 

Finally, the discussion on how the classical solution emerges from the semiclassical analysis 
is performed solving the Hamilton- Jacobi equation: one encounters a semiclassical solution 
oscillating nearby the classical solution. For A > and a ^> 1, this leads to a de Sitter 
space-time, as expected from GR. 

One then concludes that quantum cosmology applied to HL gravity suggests that this pro- 
posal matches the expectations of a quantum gravity model for the very early universe, as it 
provides, for instance, a hint for the singularity problem for the A = case. In what concerns 
specific solutions, the model suggests that GR behaviour is recovered at the semiclassical limit. 
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